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We investigate an extended version of the periodic Anderson model (the so-called periodic 
Anderson-Hubbard model) with the aim to understand the role of interaction between conduction 
electrons in the formation of the heavy-fermion and mixed- valence states. Two methods are used: (i) 
variational calculation with the Gutzwiller wave function optimizing numerically the ground-state 
energy and (ii) exact diagonalization of the Hamiltonian for short chains. The /-level occupancy 
and the renormalization factor of the quasiparticles are calculated as a function of the energy of the 
/-orbital for a wide range of the interaction parameters. The results obtained by the two methods 
are in reasonably good agreement for the periodic Anderson model. The agreement is maintained 
even when the interaction between band electrons, Ud, is taken into account, except for the half-filled 
case. This discrepancy can be explained by the difference between the physics of the one- and higher 
dimensional models. We find that this interaction shifts and widens the energy range of the bare 
/-level, where heavy-fermion behavior can be observed. For large enough Ud this range may lie even 
above the bare conduction band. The Gutzwiller method indicates a robust transition from Kondo 
insulator to Mott insulator in the half-filled model, while Ud enhances the quasi-particle mass when 
the filling is close to half filling. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.30.Mb 



I. INTRODUCTION 

Rare-earth materials exhibit numerous remarkable 
phenomena such as heavy-fermion behavior, valence fluc- 
tuations, and unconventional superconductivity. The 
simplest model that can account for these phenomena 
is the periodic Anderson model (PAM), where mobile 
conduction electrons in a broad band of width W can 
hybridize with immobile /-electrons sitting at the lattice 
sites. The Coulomb repulsion is taken into account be- 
tween the /-electrons only. Written in a mixed, Bloch 
and Wannier representation, this model is denned by the 
Hamiltonian 

k,a j.cr 
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where d ka (d k ) is the creation (annihilation) operator 
of conduction electrons with wave vector k and spin tr, 
while (fj a ) denotes the creation (annihilation) oper- 
ator of /-electrons at site rj in an arbitrary dimensional 
lattice with N lattice sites, nj a = fj C fj„ is the num- 
ber operator of /-electrons at site rj , and fij a is denned 
similarly. The hybridization matrix element between /- 
and d-states is denoted by V, and Uf is the strength of 
the on-site Coulomb repulsion between /-electrons. We 
consider the nondegenerate case, i.e., one d- and one /- 
orbital per site is assumed. Therefore, owing to the two 
possible orientations of the spin, the average number of 
d- and /-electrons per site, rid and n/, respectively, can 



vary between zero and two. The filling will refer to the 
ratio of the total electron density per site (rid + rif) and 
the maximally allowed electron density (n max = 4). 

Although it has been investigated for several decades ji 
this model and its extended versions are still in the 
forefront of condensed-matter physics. Since exact re- 
sults are available only for certain special cases r- besides 
the large number of perturbative studies nonperturbative 
techniques have also been developed to go beyond the 
weak-coupling limit. The Gutzwiller variational method- 
has been applied by several authors^ - — In this method, 
an uncontrolled approximation (the so-called Gutzwiller 
approximation^) is often used to calculate expectation 
values with the correlated wave function. Metzner and 
Vollhardtii have shown that the expectation values can 
be evaluated exactly in one dimension. Later they con- 
sidered the limit of large dimensions^ where analytic 
treatment is possible. Gcbhard- 3 developed a technique 
to calculate expectation values in a controlled expansion 
in the inverse of the degeneracy of the /-level and in the 
inverse of the dimension of the lattice. He showed that 
the Gutzwiller approximation provides exact results in 
the limit of large dimensions. Moreover, in this limit this 
method is equivalent to the slave-boson mean-field theory 
of Kotliar and Ruckenstehi i 14 ! 15 Later on, the dynamical 
mean-field theory^ which, too, is exact in the limit of in- 
finite dimensions, has been applied to the periodic Ander- 
son model by several authors^— to better understand 
the main features of the model. To avoid the problem 
related to the Gutzwiller approximation, Shiba^ applied 
the variational Monte Carlo method. The PAM was in- 
vestigated also by using the projector-based renormal- 
ization method^ for arbitrary degeneracy of the /-level. 
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The ferromagnetic properties of the PAM have been stud- 
ied with the density-matrix renormalization group32 

In view of the widespread application of the Gutzwiller 
approximation, it is important to know how reliable this 
method is. As will be demonstrated in this paper by com- 
paring the results with those of exact diagonalization, the 
Gutzwiller method gives - in spite of its limitations - re- 
liable results for the number of electrons occupying the 
/-orbital. The /-level occupancy is a significant quan- 
tity, for it has recently been proposed 2 ^ and experimen- 
tally verified^ that the pressure induced enhancement of 
the superconducting transition temperature of Ce based 
compounds, CeCu2(Ge,Si)2 is closely related to a sharp 
change of the valence of Ce. 

Several extensions of the periodic Anderson model 
have been considered so far in order to make the model 
more realistic. It was found that nearest-neighbor in- 
teraction between /-electrons affects the stability of the 
magnetic ground state in the Kondo regime^ On the 
other hand, the on-site interaction between d- and /- 
electrons (Udf Y^,j a cr' ^ja^ja') influences drastically the 
occupation number of /-electrons^ It has been shown 
that a large Udf destroys the Kondo state and narrows 
the intermediate valence regimej 23 i 26 Its treatment in the 
framework of the Gutzwiller method is, however, quite 
cumbersome. In our previous work^i we assumed a spe- 
cial form for this interaction, Udf^jh^n^n^h^, and 
pointed out that the intermediate-valence regime is nar- 
rowed in the presence of this interaction. 

The model we study in the second part of this pa- 
per includes the interaction between conduction elec- 
trons (d-electrons). Although the corresponding impu- 
rity problem has been examined thoroughly in several 
papers r^r— only a few results are available on the lat- 
tice problem.— ~— Fulde and coworkers^ have pointed 
out that the hcavy-fcrmion properties^ of Ce-dopcd 
Nd2Cu04 cannot be explained without taking correla- 
tions between conduction electrons into account. Al- 
though it has been shown^ that correlations between 
conduction electrons may increase the effective mass, and 
the competition between Coulomb repulsion in the d- 
and /-electron subsystem may lead to a transition from 
Kondo to Mott insulator, the role of the electron-electron 
interaction in the conduction electron subsystem is not 
fully clarified. In this paper we calculate the number of 
/-electrons per site and the probability of double occu- 
pancy of /-orbitals as a function of the energy of the bare 
/-level, the hybridization, and the /-/ and d-d Coulomb 
interactions. The calculations are carried out for a wide 
range of parameters of the model Hamiltonian, and the 
regions for Kondo- like behavior as well as for valence fluc- 
tuations are determined. 

The paper is divided into two main parts. Firstly, we 
investigate the reliability of the Gutzwiller method. We 
compare the variational results with those of exact diag- 
onalization on finite chains. Secondly, we analyze what 
happens when the interaction between conduction elec- 
trons, Ud, is switched on. 



II. VARIATIONAL CALCULATION AND 
EXACT DIAGONALIZATION 

A. Variational calculation 

First of all, following Ref. [36j we summarize briefly 
the main steps of the variational calculation for the origi- 
nal periodic Anderson model without interaction between 
conduction electrons, Ud = 0. In this paper we restrict 
ourselves to the paramagnetic case, i.e., the number of 
up-spin electrons, iV-f , equals that of down-spin electrons, 
N±. Furthermore, we carry out the explicite calculation 
only for the system being half-filled or less than that, 
since the results for the system more than half-filled can 
be obtained straightforwardly owing to the electron-hole 
symmetry. 

The trial wave function is chosen in the form 



i*)=^nnh/L 



v kd{ a 



|0>, 



(2) 



where the mixing amplitudes Uk and Vk are variational 
parameters. P/ is the Gutzwiller projector for /-elect- 



pf _ 



/ a/ 



(3) 



where the variational parameter r\f is controlled by Uf. 
We use the Gutzwiller approximation to evaluate the ex- 
pectation values. Optimizing with respect to the mixing 
amplitudes, we obtain 
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for the ground-state energy per site, where n/ and Vf 
denote the number of /-electrons per site and the density 
of doubly occupied /-sites, respectively, V = V \JqJ is the 
renormalized hybridization amplitude with 
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while the renormalized energy of the /-level, if has to 
be determined self-consistently from the condition 
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The k sum in Eqs. (J4j> and ((6|) [and later on, Eqs. (fT9]) and 
(pH]) in the next section] extends over the Uf = Ud = 
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Fermi sea in a manner familiar from the periodic Ander- 
son model^ since the Gutzwiller method respects Lut- 
tinger's theorem and leaves the Fermi volume unchanged. 

The quantities rif and Vf, and thereby if and qf de- 
pend on the as yet undetermined variational parameter 
rjf. Optimizing with respect to this parameter is equiva- 
lent to minimizing the energy with respect to rif and Vf, 
which leads to 
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These equations have to be solved together with the self- 
consistency condition 

The summation over k in Eqs. ^ and ([5]) could be 
carried out numerically for a realistic dispersion curve 
£d{k), but the variational procedure, i.e., the numerical 
optimization of the ground-state energy with the self- 
consistency condition (|6|), would be very cumbersome. 
Instead of that, we assume a constant density of states, 
p(e) = 1/W, in the interval e £ [-W/2, W/2], since then 
the energy density and the self-consistent value of if can 
be expressed analytically from Eqs. (|4]) and ([6]) as a func- 
tion of rif and Vf. 

However, the self-consistent solution of the minimum 
conditions for rif and Vf can be found analytically only 
in special cases, e.g., for Uf — > oo, when V <C W ^ In 
this paper we will solve Eqs. ([7]) and ((8]) numerically for 
various values of V, Uf, and Ef in order to determine 
the range of parameters for the Kondo or intermediate- 
valence behavior, and for the crossover regime between 
them. 

Firstly, we calculate the U /- and e /-dependence of the 
/-level occupancy, n/, and of the renormalization factor 
qf in the half-filled case, where the total number of elec- 
trons equals the sum of the number of d- and /-orbitals 
(the electron density per site n = rid + rif = 2), and in 
the 1/3-filled case (n = 4/3). Other fillings will be dis- 
cussed later in the next subsection, where we compare 
the results with those obtained by exact diagonalization. 

We note that our model with n electrons can be 
mapped onto a model with n holes (4 — n electrons), 
provided that the energy level of the /-hole is chosen as 
— (sf + Uf). Therefore, the results for n > 2 can be ob- 
tained straightforwardly from those for n < 2. Owing to 
this symmetry in the special, symmetric half-filled case, 
when n — 2 and the bare /-level is located at Ef = —Uf/2, 
both rif and rid are exactly equal to 1. 

The /-level occupancy is displayed as a function of 
the bare /-level energy and of Uf for V/W = 0.1 in the 
half-filled and 1/3-filled cases, respectively, in Figs.[T]and 
[2J Five different regimes can be distinguished. When 
Ef + Uf lies below the conduction band, all electrons oc- 
cupy /-orbitals, n/ « 2 and 4/3, respectively. The re- 
gions, where rif varies smoothly, almost linearly from 2 
(or 4/3) to 1 and later from 1 to 0, are the intermediate- 
valence regimes. On the plateau between them, rif devi- 



ates from unity by an exponentially small amount. This 
is, as we will see, the Kondo regime, since the double- 
occupancy rate is exponentially small here. Finally, when 
Ef lies well above the conduction band, all electrons oc- 
cupy states in the conduction band, and rif m 0. There 
are no sharp boundaries between these regimes; narrow 
crossover regions separate them. 




FIG. 1: The /-level occupancy as a function of ej and [//at 
half filling (n = 2) for V/W = 0.1. 




FIG. 2: The /-level occupancy as a function of ej and Uf at 
1/3 filling (n = 4/3) for V/W = 0.1. 

The boundary of the rif ~ 1 plateau could be defined 
by setting a somewhat arbitrary criterion for the devia- 
tion of rif from unity. Figure |3] shows nf in the Uf-Ef 
plane for a particular value of V/W using a color code. 
The "boundary" of the plateau defined by |1— n/| = 0.005 
is drawn with a white line. As can be seen in the figures, 
a plateau develops only when Uf exceeds a not sharply 
defined threshold value, Uj, which itself depends on V 
and on the total electron density. Besides V/W = 0.1, 
we have done calculations for V/W = 0.05 and 0.2, and 
obtained similar results. The upper and lower limits of Ef 
between which the plateau forms can be estimated from 
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the numerical data to be roughly 

-U } + E F {n d ) +aA f <e f < E F (n d ) - aA f , (9) 

where E F (rid) is the Fermi level of the conduction band 
with rid =n—l electrons, A/ = TipV 2 is the width of 
the /-level in the impurity problem, and a is a numerical 
factor of order 10, which depends weakly on V,Uf, and 
n. The factor a is smaller by about 10% for n = 4/3 than 
for n — 2, which shows that the plateau slightly expands 
as the filling of the conduction band decreases from half 
filling. 




£f/W £ f /W 

FIG. 3: (Color online) The /-level occupancy as a function 
of Ef and Uf at half filling (left) and 1/3 filling (right) for 
V/W = 0.1. The color code is shown at the right edge of the 
panels. The boundary of the rif ~ 1 plateau is drawn with a 
white line. 

These results are somewhat surprising. One could 
argue, based on the results for the impurity Anderson 
model that a Kondo-like behavior (i.e., « 1 with 
very small valence fluctuations) is realized when the 
Fermi level is located between the bare /-level (e/) and 
the energy Ef + Uf of a second /-electron occupying 
the same site. That is, we could expect the condition 
-Uf + E F < E f < E F , when A f < W. Condition © 
obtained by the Gutzwiller method indicates that the 
Kondo-like behavior is realized in the periodic Anderson 
model in a much narrower interval for Ef. This will be 
confirmed later by exact diagonalization. 

The /-electrons are strongly correlated on this plateau, 
since not only the average occupancy of the /-orbital is 
close to unity there, but the number of empty or doubly 
occupied /-orbitals is almost negligible. Correlations be- 
tween /-electrons can conveniently be characterized by 
the renormalization factor qf, which is simply related to 
the double-occupancy rate Vf as 

g / =8i/ / (l-2i/ / ), (10) 

when nf is exactly one. This quantity is plotted versus 
£ f and U f in Fig. H for V/W = 0.1 at half filling. 

It is clearly seen that qf decreases rapidly from about 
1, when the /-level is doubly occupied or empty, to about 
as nf approaches one from either side. When qf rj 0, 
the double-occupancy rate is also close to zero, and the /- 
electrons show heavy-fermion behavior; the effective mass 
becomes large as m* oc qj . We can, therefore, define the 
Kondo regime by setting a limit on qf, by requiring, e.g., 
qf < 0.005. This boundary is marked by a white line in 




FIG. 4: The kinetic energy renormalization factor of /-elec- 
trons as a function of Sf and Uf at half filling for V/W = 0.1. 

Fig. where qf is shown for n = 2 and n = 4/3 using a 
color code. 
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FIG. 5: (Color online) The parameter qf is displayed for 
V/W — 0.1 using a color code shown at the right edge. The 
boundary of the Kondo regime defined by qf — 0.005 is drawn 
with a white line. Left: half-filled case, right: 1/3-filled case. 

The Kondo regime thus defined appears again above 
a critical Uj, which is, however, somewhat larger than 
the one found earlier, since the criterion |1 — n/| < 0.005 
is less strict than the condition qf < 0.005. In this lat- 
ter case the probability of double occupancy has to be 
less than 0.0006. Nevertheless, comparison with Fig. [3] 
shows that apart from a rounding around the critical Uf, 
the two criteria define the same regime. The plateau 
slightly expands when the filling of the conduction band 
decreases. 

When nf is exactly one, and Eq. (fTU)) holds, Eq. (0 
can be easily solved in the limit V <C W. We get 

nd ( u f \ /-i-n 

qf = WIW exp V^yw)' 

where rid is the number of the conduction electrons per 
site. The factor rid m the prcfactor explains why the 
critical Uf gets smaller as the filling decreases. 

The total energy density takes a simple form in this 
limit, 

W { U f ] 

£ = ef + E d (n d ) - n d — exp |- mv y w J, (12) 
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where the first term is the energy of the half-filled /- 
orbital without polarity fluctuations, E^n^) is the en- 
ergy of the decoupled conduction band for filling rid, and 
the last term describes the coupling between /-electrons 
and conduction electrons, in other words, the energy de- 
crease owing to the polarity fluctuations caused by d-f 
hybridization. This term arising from the Kondo effect 
gives the characteristic energy scale in the Kondo regime. 
The Kondo energy, Ek , is defined by the energy decrease 
per conduction electron, i.e., 

Ek = ? < 13 > 

In the remaining part of this subsection, we study more 
quantitatively the dependence of the threshold value of 
Uf on V in the half-filled case. Figure [6] shows q/ as 
a function of V for several values of Uf at £/ = —Uf/2, 
where n/ is exactly one. The threshold values determined 

1 j i i i i i i i i i i i i i i i i i i i 




V/W 

FIG. 6: qj vs. V in the symmetric half-filled case, ej = 
— Uf/2, for Uf/W = 1 (continuous line), 3 (dashed), 5 (dot- 
ted), 7 (dashed with one dot), and 20 (dashed with two dots), 
respectively. 

from qf(Uj, V) = 0.005 are given in Table U together the 
corresponding Kondo coupling J — 8V 2 /Uj, since, in 
the Kondo regime, the periodic Anderson model can be 
mapped onto a Kondo lattice model (KLM). 



V/W 


uyw 


J/W 


0.16 


3 


0.066 


0.21 


5 


0.071 


0.26 


7 


0.076 


0.48 


20 


0.093 


0.74 


40 


0.110 


0.97 


60 


0.125 


1.17 


80 


0.136 


1.36 


100 


0.148 



TABLE I: The critical Uf of the Kondo plateau for several 
values of V and the corresponding Kondo coupling. 



The dependence of Uj on V can be fitted by the ana- 
lytic functional form 

U c f /W = &2.bQ{V/W) a . (14) 

with a = 1.54 

Since by definition there are no doubly occupied or va- 
cant /-orbitals in a KLM, a rigorous mapping from PAM 
to KLM should be possible in the limit Vf — > 0. Set- 
ting a smaller limit for qf in the criterion for the Kondo 
regime, larger exponents, given in Table [TTl and larger 
numerical prcfactors are obtained in Eq. (|14[) . The expo- 
nent a seems to converge to 2 in the limit qf — > 0, which 
means that Uj is proportional to Ay, and the propor- 
tionality factor is of order hundred instead of the factor 
a rj 10 in Eq. ©. This difference is due to the stricter 
condition on q and to the rounding of the boundary at 
the critical U% 
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a. 


1.70 


1.80 


1.83 


1.86 


1.89 


1.91 



TABLE II: The exponent a in Eq. (Ti"4)) calculated for several 
threshold value of qf. 

Sinjukow and Nolting^ have shown that in the ex- 
tended Kondo limit, when Uf — > oo and V — > oo with 
V 2 /Uf remaining finite, the symmetric periodic Ander- 
son model can be mapped exactly to the Kondo lattice 
model with finite Kondo coupling. The results obtained 
by the Gutzwiller method are in agreement with this. 

B. Comparison with exact diagonalization 

With the aim to compare the variational results with 
those of a completely different method, we also performed 
exact diagonalization on relatively short chains. In order 
to check whether the results obtained for these chains 
are representative for bulk materials, we calculated the 
/-level occupancy, nj, and the density of doubly occu- 
pied /-sites, Vf, in the nonmagnetic (Sl ot = 0) ground 
state for chains of four, five, and six sites. It turned out 
that the results were in excellent agreement with each 
other, which suggests that the six-site chain behaves al- 
most like the bulk in this respect. This is in agreement 
with the finding of Chen and Callaway^ who compared 
the ground-state energy obtained from exact diagonal- 
ization of a four-site chain with Monte Carlo result on a 
sixteen-site chain. In what follows we present the results 
obtained for a six-site chain with 12, 10, 8, and 6 elec- 
trons. The case with 6 electrons is not interesting from 
the point of view of Kondo physics, because the conduc- 
tion band is exhausted when rif = 1. Nevertheless, it is 
used in the comparison of the two methods. 

The kinetic energy of conduction electrons moving 
along the chain is described by hopping between nearest- 
neighbor d-orbitals with hopping rate t, thus the band 
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width is now At. Therefore, we identify W with At, when 
comparison with the results of the variational calculation 
is made. 




FIG. 7: (Color online) (a) The /-level occupancy vs. e/ at 
2V/W = 0.375. The curves are obtained by the Gutzwiller 
method, while the symbols indicate the results of exact diag- 
onalization for 2Uf/W = (black continuous line, .), 3 (red 
dashed line, ■), 5 (green dotted line, ♦), 7 (blue dashed line 
with one dot, i) and 10 (purple dashed line with two dots, t). 
(b) The renormalization factor qf. The notation is the same 
as in panel (a). 



The /-level occupancy obtained by the two methods 
are directly compared in Fig. [Zta). As for z/y, we compare 
the results indirectly, through qf. Although this quantity 
is specific to the Gutzwiller method, it shows the strength 
of correlations more visibly than Vf itself, therefore, we 
define qf with the help of Eq. ([5]) from nf and Vf obtained 
from the exact ground-state wave function. Comparison 
with the result of the variational calculation is shown in 
Fig.Etb). 

As is seen in Fig.[7]Ja), the two methods give very sim- 
ilar results as far as the "global behavior" of the /-level 
occupancy and the extent of the n/ w 1 plateau is con- 
cerned, even though the density of states is not identical 
in the two calculations. This indicates that Eq. (|9|) found 
in the Gutzwiller method for the boundary of the Kondo 
regime is not due to the Gutzwiller approximation, but is 
a consequence of strong correlations in the lattice model. 

We find a subtle difference, however, in Fig. [TJb), 



where q/ is plotted as a function of Sf. One sees that 
qf calculated in the Gutzwiller method approaches zero 
faster in the Kondo regime than that provided by ex- 
act diagonalization. The former exhibits the exponential 
behavior given in Eqs. (fTTj) and (fT3f typical for Kondo 
physics, while the latter cannot be fitted to such a curve. 
We will discuss this quantitatively later on. 

Next we check the dependence of the Kondo plateau 
on the strength of the hybridization. In Fig. Ufa) we plot 
nf as a function of £/ for three values of V/W in the 
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FIG. 8: (Color online) (a) The /-level occupancy vs. £/ at 
Uf/W = 2.5. The curves are obtained by the Gutzwiller 
method, while symbols denote the results of exact diagonal- 
ization for 2V/W = 0.2 (black continuous, .), 0.375 (blue 
dashed, ♦), and 0.7 (red dotted, ■), respectively, (b) The 
renormalization factor qf. The notation is the same as in 
panel (a), (c) n f vs. V at Uf/W = 2.5, for e f /W = -0.75 
(black continuous line, .) and ef/W = —0.25 (blue dashed 
line, ■). 
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V/W v/W 

FIG. 9: The double-occupancy rate of /-electrons vs. V/W . 
The dotted curves indicate the results of exact diagonal- 
ization, while the continuous curves are calculated by the 
Gutzwiller method. Uf jW is 5 and 100 in panel (a) and 
(b), respectively. 



half- filled case. It is clearly seen that the plateau (i.e., 
the Kondo regime) rapidly shrinks as V increases, and 
disappears, in agreement with the results presented in 
the previous subsection. Figure [5Jb) , where qj is plot- 
ted, shows directly the disappearance of heavy-fermion 
behavior. Finally nj is plotted MS i\ function of V in 
Fig. [SJc) for two values of Sf/W. We find again that 
the two methods yield similar results for n/, but the In- 
dependence is different near the boundary of the Kondo 
regime. 

In order to better see this difference, we calculate the 
double-occupancy rate of /-electrons in the symmetric 
{sj = —Uf/2) half-filled case as a function of V near the 
boundary of the Kondo regime, i.e., where Vf <C 1. We 
find, as seen in Fig. |H1 that in contrast to the results of 
the Gutzwiller method, the dependence of Vf on V 2 /Uf 
is not exponential; Vf varies as a power of V 2 /Uf. 
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FIG. 10: (Color online) (a) The /-level occupancy vs. ej 
at Uf/W — 2 for different fillings. The hybridization is 
V/W = 0.1 in all cases. The curves are obtained by the 
Gutzwiller method, while the symbols are the results of ex- 
act diagonalization. The number of electrons per site is n — 2 
(black continuous line, .), 5/3 (red dashed line, ■), 4/3 (green 
dotted line, ♦), 1 (blue dashed line with one dot, a), respec- 
tively, (b) The /-level kinetic energy renormalization factor. 
The notation is the same as in panel (a). 



W ( V 2 \ W ( V 2 \ 2 

»> = A v l {wv,) +B v l {wv}) ■ (15) 

where A is close to unity and B w 50. This power-law- 
like dependence may be due to the small system size in 
the exact diagonalization. 

Finally we study the filling dependence of the Kondo 
regime. The /-level occupancy is shown for several fill- 
ings in Fig. [Tola). The overall agreement between the 
two methods persists as we move away from half filling, 
though its degree varies somewhat, e.g., the agreement 
in the n = 2 or 5/3 case is noticeably less good than for 
n = 4/3. This indicates that the Gutzwiller type para- 
magnetic wave function is more appropriate for metallic 
systems with few conduction electrons than for insula- 
tors. The Kondo plateau shifts towards lower /-level en- 
ergies as the filling decreases owing to the decrease of the 
Fermi level. We find a similar slight difference between 
the results of the two methods displayed in Fig. [TUJb), 
where qf is plotted as a function of £/. 



III. THE ROLE OF INTERACTION BETWEEN 
CONDUCTION ELECTRONS 

A. Variational calculation 

As a next step, we consider what happens when the 
interaction between conduction electrons is switched on. 
For the sake of simplicity a local, on-site interaction is 
assumed and the Hamiltonian takes the form 

n = Hpam + u d J2 > ( 16 ) 



where 'Hpam is the PAM Hamiltonian defined in Eq. ([T]) 
and Ud is the strength of the Coulomb interaction be- 
tween conduction electrons. This model is also known 
as the periodic Anderson-Hubbard model. At half filling 
the symmetric model corresponds to e/ = —Uf/2 + Ud/2, 
where nt — rid = 1- 

The variational calculation can be performed by a sim- 
ple generalization of the procedure used for Uf — ^ ooi2s 
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The trial wave function is chosen in the form 



I*' 



= %Pg 



10), (17) 



k a 



where Pq contains the variational parameter rjf, and an 
extra Gutzwiller projector has been introduced for d- 
electrons, which is written as 



(18) 



The variational parameter r\ d depends on E/j. Performing 
the optimization with respect to the mixing amplitudes 
we get 



fceFS 



£ = Jf H qd£d ^ + ~ e i ~ \/[<lded(k)-E f ] 2 +W' 



+(ef - gf)nf + U d v d + Ufv f 



(19) 



for the ground-state energy density, where v d is the den- 
sity of doubly occupied <i-sites, and q d denotes the kinetic 
energy renormalization factor of d-electrons given by 



1 



Id 



\ 1 _ «J| "d 
V 1 2)2 



Ud 

2 



Vd Vd 



+ ]J (y - Vd) (1 - n d + v d ) 



(20) 



which is formally identical to that found in the Hub- 
bard models The renormalized hybridization amplitude 
is now V = V \jq d qf] the other notations are the same as 
in the previous section, and the self-consistency condition 
[see Eq. ©] is now given by 



N ^ 



fceFS 



qd£d(k) - if 



[q d e d (k) -if] 2 +AV 



(21) 



The summation over k and the numerical optimization 
of the energy density with respect to nf, i/f, and 1/4 are 
carried out in the same way as in the previous section. 
The equations determining nf,Vf and Vd are now 



-u d 



d£ dqd d£ dqf 

dqd dnf dqf dtif 

d£ dqf d£ dif 

dqf dvf dif duf' 



d£ di 



dif dnf 



d£ dq d 
dqd dv d 



d£ dif 
di f dv d ' 



, (22) 



(23) 



(24) 



First we derive analytic results from these equations in 
the weak hybridization limit up to 0((V/W) 2 ) for arbi- 
trary Uf at special fillings: for nf = 1 and n d arbitrary; 
for nd = 1 and rif arbitrary; and finally for nf = rid = 1. 
Similar results were obtained in Ref. [3r| . but only for 
Uf -»• 00. 



We know that the interaction between conduction elec- 
trons suppresses charge fluctuations in the Hubbard sub- 
system. This influences the Kondo physics in the follow- 
ing ways: 

(i) Ud shifts the Fermi energy of the conduction band. 
For nf = 1 and nd < 1 we get 



E ¥ (n dl U d ) 
+ W 



rid 1 1 
T - 2 ] QdW 



Ud 

2 



W) q d 



F 1 - (25) 
dn d 



The third term in the square brackets is the contribution 
of d-f hybridization. Without it we recover the equation 
determining the Fermi energy of the Hubbard model for 
filling nd- Note that the values of Vf and v d in qf and q dl 
respectively, should be taken from the solution of Eqs. 
([2U| and ([M]) . Equation (|2"5j) has no simple closed form 
solution for arbitrary Ud except for the half-filled case, 
where E-p{nd = 1, C/d) = Ud/2- At other fillings we can 
expand E-pfrid, Ud) in the weak- or strong-coupling limit 
(U d < W or U d > W) as 

E F {n d , U d ) PS E F (n dl 0) + ^-U d + 0{U 2 d /W) (26) 

or 

E F (n d , U d ) w E F (n d , 0) + y W + 0(W 2 /U d ), (27) 

respectively. The shift of the Fermi energy is at most 
W/6 for n — 4/3 (i.e., n d = 1/3), which is much smaller 
than the shift in the half-filled case. 

(ii) Switching on Ud reduces the Kondo energy^ When 
nf = 1, we can calculate qf and the total energy density 
for finite Ud and for arbitrary n d (assuming V W). 
Instead of Eqs. (JTT]l and {12]) we find 



Qf 



n d q d 



A{V/Wf 



cxp 



U f 



wv 2 /w 



(28) 



and 



£ = Ef + E d (n d , U d ) 



W 

ndqd-^- exp 



16V 2 )w\ (29) 



where the second term of the right hand side is the energy 
of the decoupled correlated conduction band. Compared 
with Eqs. (fTTj) and (|T3"|) . qf and the exponential Kondo 
scale are reduced by q d , which is rather small when n d w 
1 (see below). For n slightly less than the half-filled case 
this mechanism yields a significant mass enhancement. 
We get q d ~ 1/5 for n = 1.95 and U d = 2 AW, which 
means that the effective mass is five times bigger for these 
parameters than without Ud- 

(iii) The most interesting effect of Ud is the Mott tran- 
sition which occurs in the Hubbard model at half filling 
(n d = 1). In the Gutzwiller-type treatment of U d it is 
known as the Brinkman-Ricc transition. It occurs when 
qd becomes zero for a finite Ud- A similar transition may 



9 



take place in the half-filled periodic Anderson-Hubbard 
model. In this model, however, even if n = 2, the Kondo 
physics may compete with Mott physics, rid and rif de- 
pend on Ud, Uf, V, and £f owing to the d-f hybridiza- 
tion, and the conditions for the Mott transition may not 
be so simple as in the Hubbard model. In what follows we 
first show in the framework of the Gutzwillcr treatment 
that the necessary conditions for the Mott transition is 
that both the /- and d-electron subsystem be half filled, 
i.e., nd = 1 and rif = 1 be fulfilled simultaneously, and 
moreover the system be in the Kondo regime. 

We see from Eq. (|20p that qd is zero only when rid = 1 
and Vd = 0. Similarly it follows from Eq. ([5]) that <?/ 
vanishes only if rif = 1 and Vf = 0. When rid = 1, the 
renormalization factor qd is simply 8^d(l — 2v d ), and Eq. 
(|M]) gives 



Ud 
W 



4 \WJ q d 



8(1 - Av d ) = (30) 



for V <C W and rif arbitrary. The second term in the 
square brackets is the contribution of d-f hybridization. 
Without it we recover the equation determining the op- 
timum Vd of the half-filled Hubbard model. 

It follows from this equation that Vd goes to zero as 
Vd approaches a finite critical value only if qf also ap- 
proaches zero, and qf is of the same order as qd- This 
situation can be realized only if rid = 1 and n/ = 1 arc 
simultaneously fulfilled, and moreover qf is given by Eq. 
(|28|) . i.e., the system is in the Kondo regime. 

When rid = 1 and rif = 1 are simultaneously satisfied, 
and the system is in the Kondo regime, the term in (|30[) 
due to d-f hybridization is independent of Vd and is equal 
to E K /W [see Eqs. (HI, (EHJ), and (|29])]. Equation ([30]) 
is easily solved to give 



V d = T 



U d 



8(W + AE K ) ' 



(31) 



which shows that Vd decreases linearly as Ud increases 
and reaches zero at U% = 2(W + AEk)- At this value 
of Ud, which - owing to the coupling between the d- and 
/-electron subsystems - is slightly larger than the critical 
value in the Hubbard model (U% = 2W), the conduction 
band undergoes a Brinkman-Rice transition. Note that 
the exponentially small correction has been neglected in 
Ref. [361 ] . Since Vd = Vf = at this transition, all po- 
larity fluctuations are suppressed, the effective d-f hy- 
bridization (V = V^/qdqj) as well as the Kondo energy 
scale become zero, that is, the Kondo effect is completely 
quenched. The system transforms from a Kondo insula- 
tor into a Mott insulator. 

Analytically we can claim only that the condition rid = 
rif = 1 is realized in the symmetric point of the half- 
filled model, where £/ = —Uf/2 + Ud/2. Indeed, when 
Ud is smaller than U% and is not very close to it, it is 
found numerically that rid = nf = 1 is realized only at 
the symmetric point, and thus one could expect that a 
Brinkman-Rice transition occurs only in the half-filled 



symmetric periodic Anderson-Hubbard model, and that 
the system becomes a Mott insulator for Ud > U\ only if 
e f = -U f /2 + U d /2. 

Contrary to this expectation we have found numeri- 
cally that when Ud is slightly smaller than the critical 
value, rif = rid = I holds not only at the symmetric 
point, but - within the limits of the numerical accuracy of 
our calculations, which was about 10 -6 - in a wide range 
of £f within the rif « 1 plateau. In order to find the ex- 
tent of this range, we display the Ud- and £/-dependcnce 
of the /-level occupancy and of the renormalization fac- 
tor q f for U d < 2W at half filling in Figs. [TT] and [T2] 
respectively. 




FIG. 11: The /-level occupancy vs. the /-level energy and Ud 
at half filling for V/W = 0.1 and U f /W = 2. 




FIG. 12: The kinetic energy renormalization factor for /- 
electrons vs. ej and Ud at half filling for V/W = 0.1 and 
Uf/W = 2. 

It is clearly seen that the Kondo plateau, where rif rj 
1 and qf « 0, shifts towards higher energies owing to 
the shift of the Fermi energy by Ud/2, and the center of 
the plateau is located indeed at Sf = — Uf/2 + Ud/2 as 
expected from Eq. (j25]) . The condition for the Kondo 
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regime can be written similarly to Eq. J9]) as 

-Uf+E F (n d ,U d )+aA f <£/ < E F (n d ,U d )-aA f . (32) 

Note that the center of the plateau is at the center of the 
nonintcracting d-band, when U d = Uf. In other words, 
the /-level does not need to lie low enough compared to 
the conduction band to show heavy-fermion behavior. 

Another remarkable feature is that the plateau widens 
as U d approaches the critical value U^. At U d = 2W, it 
is situated in the range — Uf + U d /2<Sf<U d /2. That 
means that the narrowing of the plateau compared to the 
impurity model given by a A f in Eq. Q gets remarkably 
smaller close to U d . This is probably due to the formation 
of the Hubbard subbands and the drastic variation of the 
density of states at the Fermi energy near the transition 
point. 

Numerical calculations give n d = Uf = 1 at U d on 
the whole Kondo plateau. This indicates that - at least 
within the Gutzwillcr-type treatment of correlations - 
both n d and rif are fixed to exactly unity in the half- 
filled model as we approach U d and the condition for 
Kondo behavior is satisfied. The renormalization factors, 
q d and qf vanish simultaneously, the d-f hybridization is 
completely suppressed, so is the Kondo effect, and a Mott 
transition takes place. This transition in the conduction 
electron subsystem is robust, it is the dominant feature 
of the half-filled model. 

Our finding that the Brinkman-Ricc transition and the 
Mott insulating state are not restricted to the symmetric 
model is corroborated by calculations for U d > U d . The 
numerical variational calculation yields meaningless neg- 
ative values for v d in the whole interval —Uf + U d /2 < 
Ef < U d /2. Note that for £f outside this interval we can 
carry out the numerical calculations for arbitrary large 
U d without any difficulty. 

Next we show that the d-f hybridization prevents the 
Mott transition when n < 2 (or for n > 2). In this case 
the term coming from the d-f hybridization in Eq. (|30[) 
becomes large, if q d — > 0, since qf is always finite for 
Uf < 1, and thus there exists no such solution for v d (or 
q d \ which approaches zero at a finite U d . Charge fluc- 
tuations on the ci-orbitals are thus not completely sup- 
pressed. A finite v d indicates the existence of a Fermi 
surface, since q d is identified with the discontinuity at 
the Fermi wavenumber in the single-particle occupation 
number^ This can be understood as follows: even if the 
correlated conduction band is half filled and U d is large 
enough, so that the conduction band is separated into 
Hubbard subbands and the Fermi level lies within the 
/-band located in the Hubbard gap, the d-electrons are 
taking part in the formation of the Fermi surface via d-f 
hybridization. 

The results of the numerical calculations in the 1/3- 
filled case (n = 4/3) are shown for < U d < 3W in Fig. 
[T3l We observe that one more plateau appears at higher 
/-level energies, above the bare conduction band, when 
U d > 2W, corresponding to n d ^ 1. Its formation in- 
dicates that two separate Hubbard subbands arc formed 



above this critical value of U d . The plateau appears when 
Ef is located between the two subbands, i.e., in the Hub- 
bard gap. The Fermi level is located in the /-band in this 
situation. The center of the plateau is approximately at 
U d /2, which indicates that the upper and lower subbands 
are centered at and U d , respectively, i.e., the location 
of the subbands is the same as in the Hubbard model. 




FIG. 13: The /-level occupancy vs. /-level energy and Ud at 
1/3 filling (n = 4/3), for V/W = 0.1 and Uf/W = 2. 

We observe furthermore that the n/ « 1 plateau hardly 
shifts as U d increases. Since according to Eq. (f2~T| the 
Fermi energy only weakly depends on U d away from half 
filling, we find again the condition given in (|32[) for the 
Kondo regime. 

It is worth mentioning here what happens when the 
system is more than half filled. The answer can be ob- 
tained without any further calculation from electron-hole 
symmetry. A model with n electrons can be mapped onto 
a model with n holes (4 — n electrons) by the transfor- 
mation: 

dl^e^d hjs , d ja ^e-^dt j9 , 
fl^-e^hft, fir^-e-*'^, (33) 

where the index (h) refers to holes, and a = — a. If 
the kinetic energy of conduction electrons is written in 
Wannier representation, 

E^( fe )4<A CT = E^A' ( 34 ) 

k.a ijo 

and the phase factor is chosen in the form ipj = Q ■ Tj , it 
is easily seen that the kinetic energy term is transformed 
into 

J2 £hd( - k ^lkJhk*> ( 35 ) 

where 

£hd(fc) - -e d {k + Q). (36) 
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Assuming that ta = 0, the center of the band sets the 
zero of energy. The term describing hybridization is in- 
variant under this transformation, while the on-site en- 
ergy of /-levels and the on-site interaction terms give rise 
to energy shifts. Therefore the Hamiltonian written in 
terms of the creation and annihilation operators of holes 
has the same form as for electrons with shifted energies 
for the d- and /-electrons, and an overall energy shift: 



n c (e d (k),e f ,V,U d ,U d ) — > 
H h (e hd (k) - U d , e hf ,V, U d , U f ) + E , 



(37) 



where e hf = -e f -U f and E = (2e } + U } + U d )N . If the 
energy levels are measured from —U d the Hamiltonian in 
hole representation becomes 

H h (e hd (k), ehf , V, U d , Uf) + E , (38) 

with e"h/ = —£f — Uf + U d and 

E = -(2e d + U d )N h + (2e f + U f + 2e d + U d )N, (39) 

where A h is the total number of holes. Provided that 
Ei ld (k) = —e d (k + Q) = e d {k) for a certain Q, as is the 
case for the one-dimensional model with nearest-neighbor 
hopping, or when a constant density of states is assumed, 
then the dispersion curve of d-holes is the same as for d- 
electrons and the results obtained in the electron repre- 
sentation can be applied to holes when the energy shifts 
are taken into account. 

Using this transformation, the results for n > 2 can be 
obtained straightforwardly from those for = 4 — n < 
2. We can get, for example, the Fermi energy of the 
correlated conduction band for n d = n — 1 > 1 (n/ = 
1) from that for n^ d = 2 — n d < 1 [see Eq. ((25])] by 
first shifting the origin of the energy by —U d , and then 
reversing the energy axis. We get 

E F (n d ,U d ) = -[E F (n hd ,U d )-U d ], (40) 

from which for n d > 1 

E F (n d ,U d = 0) = -E F (2-n d ,U d = 0). (41) 

The equation giving the shift of the Fermi energy for 
n d > 1 is thus 

E F (n d , U d ) « E F (n d , 0) - + U d + 0(W 2 /U d ) 

(42) 

instead of Eq. ([2T]) . This shows that the shift of the 
Fermi energy owing to U d for n d > 1 is larger than that 
at half- filling. 

The condition on for the Kondo regime is obtained 
for n > 2 as follows: The condition on the /-hole level 
£h/ for ??h < 2 is formally the same as for electrons [see 
Eq. ([32]) ]. since the Hamiltonian has the same form, i.e., 



The condition on the /-electron level for the Kondo 
regime for n > 2 is simply obtained by rewriting this 
condition for the original Ef using £hf = — £/ ~ Uf + U d . 
We get 



-U f + U d - E F (n hd , U d ) +aA f < 
£ f ^U d - E F (n hd , U d ) 



aA 



(44) 



U f + E F (n hd , U d ) + aA f < e hf < E F (n hd , U d ) - aA f . 

(43) 



Since according to Eq. (|40]) U d — E F (rihd,U d ) (n\ id < 1) 
is the Fermi energy of the interacting conduction band, 
E F (n d , U d ), for n d = 2 — rih d > 1, the condition takes the 
same form given in Eq. ([3^]) for all fillings. For n > 2, the 
shift of the n/ ~ 1 plateau is thus even larger than in the 
half-filled case, it may appear above the bare conduction 
band. 



B. Comparison with the results by exact 
diagonalization 

Now, we compare the results of exact diagonalization 
with those obtained by variational calculation. Note that 
we will discuss only the case n > 1, i.e., more than six 
electrons on a six-site chain. The quarter-filled case is 
not interesting from the point of view of Kondo physics, 
because the conduction band is exhausted when n/ = 1 . 

The values of nf obtained by both methods are shown 
for several fillings at Uf = U d = 2W in Fig. H3Ta). The 
overall agreement between the results of the two methods 
demonstrated earlier remains good for finite U d . It is 
remarkable that the agreement is even better than for 
U d = 0. 

The shift of the nf rj 1 plateau due to U d is observed 
in both methods, showing that the shift of the plateau is 
not an artefact of the Gutzwillcr approximation, and may 
be observable in some materials, where the conduction 
electrons are strongly correlated, i.e., they may exhibit 
heavy-fermion behavior despite the fact that the bare /- 
level does not lie below the conduction band. 

The formation of two separate plateaus corresponding 
to nf w 1 and n d pa 1 is also observed in both methods. 
The formation of the Hubbard subbands owing to U d in 
the conduction-electron subsystem is thus confirmed by 
exact diagonalization, too. We carried out the compari- 
son also at Uf = U d = 5W for n = 5/3, and found that 
the agreement between the two methods is almost perfect 
for n/ . 

Figurc[l5]shows the kinetic energy renormalization fac- 
tor of conduction electrons (q d ) for three different fillings. 
The agreement between the two methods arc fairly good 
in panel (a), while a marked difference is seen in panel 
(b), i.e., at half filling. At this point we should note that 
the two methods are complementary. The Gutzwillcr 
method is exact for large dimensions, while the exact 
diagonalization was performed on chains. We argue that 
the difference between the results obtained by the two 
methods at half filling is due to the unusual behavior of 
the onc-dimcnsional half-filled Hubbard model. 
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FIG. 14: (Color online) (a) The /-level occupancy vs. ej 
at Uf/W = U d /W = 2 and V/W = 0.1. The curves are 
obtained using the Gutzwiller method, while the symbols are 
the results of exact diagonalization. The electron number 
per site is 2 (half filling) (black continuous line, .), 5/3 (red 
dashed line, ■), and 4/3 (green dotted line, ♦), respectively, 
(b) The renormalization factor qj . The notation is the same 
as in panel (a). 

The one-dimensional Hubbard model can be solved ex- 
actly by Bethe ansatz4^ At half filling its ground state 
is conducting only for Ud = and insulating for any 
nonzero Ud- On the other hand, in higher dimensions the 
half-filled Hubbard model is expected to remain metallic 
until a finite critical where the Mott transition takes 
place. Therefore, when we compare the results obtained 
by the Gutzwiller method and by exact diagonalization 
of the Hamiltonian of a chain displayed in Fig. [PsT b) 
and interpret the difference, we have to keep in mind the 
fundamental difference between the physics of one- and 
higher dimensional half-filled Hubbard models. 

The Gutzwiller method gives not only a vanishing va- 
lence fluctuation on /-orbitals, vt = 0, and consequently 
qj- = at Ud = 2W + 8Ek, but the same is true for the 
conduction-electron subsystem: also v& and qd vanish at 
the Brinkman-Ricc transition. On the other hand, exact 
diagonalization gives a finite qd in agreement with the 
known behavior of the one-dimensional half-filled Hub- 
bard model, where Vd is finite for arbitrary Ud^ We 
believe that the disagreement between the predictions 
of the Gutzwiller method and of exact diagonalization 
seen in Fig. ITSl b) is thus the consequence of the different 
behavior of the one-dimensional and higher dimensional 
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FIG. 15: (Color online) The kinetic energy renormalization 
factors of the conduction electrons vs. e/ at Uf/W = Ud/W = 
2. The curves are obtained by the Gutzwiller method, the 
symbols show the values calculated with exact diagonaliza- 
tion using Eq. (|20[1 . (a) The results for n — 5/3 (the red 
dashed line, ■) and n=4/3 (the green dotted line, ♦). The 
hybridization is V/W = 0.1 in all cases, (b) The results for 
half filling, n — 2 (the black continuous line, .). 



periodic Anderson-Hubbard models. 



IV. CONCLUSIONS 

In this paper we considered an extended periodic An- 
derson modell, the so-called periodic Anderson-Hubbard 
model with on-site Coulomb repulsion in the conduction- 
electron subsystem. Our main aim was to investigate 
how the additional repulsive interaction between con- 
duction electrons influences the Kondo regime, and how 
the Kondo physics and Mott physics compete. For this 
study we calculated the average number of /- and d- 
electrons per site, nj and rid, and the probability of dou- 
ble occupancy in both subsystems, vj and Ud, using the 
Gutzwiller variational method. In order to check the re- 
liability of this method, we also performed exact diago- 
nalization on relatively short chains. Since to our best 
knowledge no such comparison was presented even for the 
original periodic Anderson model, we also present results 
for the original periodic Anderson model. 

A rather good agreement was found between the re- 
sults of the two methods in the original model as far as 
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the location of the Kondo and valence- fluctuation regimes 
are concerned. A subtle difference was, however, found 
near the boundary of the Kondo plateau. Namely, while 
the results of the Gutzwiller method exhibit an exponen- 
tial dependence of the double occupancy on the charac- 
teristic combination of the couplings, V 2 /Uf, those of 
exact diagonalization show a power-law behavior. This 
will be the subject of subsequent studies. 

The situation is somewhat different for the extended 
model. Both methods indicate that when the on-site 
Coulomb repulsion between conduction electrons (Ud) 
is switched on, in the half-filled case the heavy-fermion 
regime shifts towards higher energies of the bare /-level 
by Ud/2 in accordance with the shift of the Fermi energy 
owing to Ud- A marked difference appears, however, be- 
tween the results provided by the two methods, when Ud 
is of the order of 2W. The Gutzwiller method indicates 
that both Ud and rif are fixed to unity for a wide range 
of the /-level energies in the half-filled model at a crit- 
ical value of Ud, and a robust Brinkman-Rice transition 
takes place to a Mott insulator. Although the Kondo ef- 
fect is enhanced, when the Anderson-Hubbard system ap- 
proaches the critical point, this effect is completely sup- 
pressed right at the transition, and all charge fluctuations 
are suppressed. Even though the exact diagonalization 
on chains does not reproduce this result, we believe that 
the Gutzwiller method describes correctly the scenario 



in higher dimensional systems, and the different behav- 
ior found for linear chains is simply a consequence of the 
anomalies of low-dimensional systems. 

When the electron system is less than half filled, the 
Mott transition is suppressed by the d-f hybridization, 
and besides the Kondo plateau (rif sa 1) another plateau 
appears at rid ~ 1, provided that Ud > 2W, i.e., when 
the conduction band is split into a lower and upper Hub- 
bard band. The results provided by the two methods 
are in surprisingly good agreement in this case, in par- 
ticular when correlations are strong. The shift of the 
heavy-fermion regime towards higher bare /-level ener- 
gies owing to Ud is small compared to that in the half- 
filled case, because the shift of the Fermi energy due to 
Ud is at most ridW/2 for < 1. On the other hand, 
when the electron system is more than half filled, the 
shift of the Kondo regime with Ud is much larger, since 
the shift of the Fermi energy is also larger than that in 
the half-filled case. 
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